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Abstract 

We will provide a review of some of the physics which can be addressed by 
O ■ studying fluctuations and correlations in heavy ion collisions. We will discuss Lattice 

^ . QCD results on fluctuations and correlations and will put them into context with 

observables which have been measured in heavy-ion collisions. Special attention will 
be given to the QCD critical point and the first order co-existence region, and we 
^ . will discuss how the measurement of fluctuations and correlations can help in an 

! experimental search for non-trivial structures in the QCD phase diagram. 

in 

1 Introduction 

^ ■ Fluctuations and correlations are important characteristics of any physical system. They 
provide essential information about the effective degrees of freedom and their possible 
quasi-particle nature. In addition, the susceptibilities, which characterize the correlations 
^ ' and fluctuations, determine the response of the system to small external forces. 

• In general, one can distinguish between several classes of fluctuations. On the most 

fundamental level there are quantum fluctuations, which arise if the specific observable 
does not commute with the Hamiltonian of the system under consideration. These fluc- 
tuations probably play less a role for the physics of heavy-ion collisions. Second, there 
are "dynamical" fluctuations and correlations reflecting the underlying dynamics of the 
system. Examples are density fluctuations, which are controlled by the compressibility of 
the system. Finally, there are "trivial" fluctuations induced by the measurement process 
itself, such as finite number statistics, etc. These need to be understood, controlled and 
subtracted in order to access the dynamical fluctuations which tell as about the properties 
of the system. 

A prominent example where the measurement of correlations has lead to a scientific 
breakthrough are the fluctuations of the cosmic microwave background radiation [Tj, first 
carried out by the COBE satelite [2] and later refined by WMAP In the case of the 
cosmic microwave background, the fluctuations are on the level of 10"'* with respect to 
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the thermal distribution. In addition, a large dipole correlation, due to the motion of 
the earth through the heatbath of the microwave background, has been seen and needed 
to be subtracted before the more interesting correlations due to the big bang, could be 
revealed. 

In the cosmic microwave background one observes spatial correlations. In contrast, 
in heavy-ion collisions we are restricted to correlations in momenta and the connection 
to spatial correlations is, in general, non-trivial. Nonetheless, in the case of heavy-ion 
collisions we have a quite similar situation. To first order we observe thermal spectra 
and radial flow. Next, we see a large quadrupole correlation due to elliptic flow. The 
interesting question then remains, if we will be able to identify smaller correlations due to 
a possible phase transition and the QCD critical point after subtracting the two dominant 
backgrounds, thermal emission and elliptic flow. 

Fluctuations are closely related to phase transitions. The well known phenomenon 
of critical opalescence is a result of fluctuations at all length scales due to a second 
order phase transition. First order transitions, on the other hand, give rise to bubble 
formation, i.e. density fluctuations at the extreme. Considering the richness of the QCD 
phase-diagram, as sketched in Fig{Tl the study of fluctuations in heavy-ion physics should 



lead to a rich set of phenomena and is an essential tool for the experimental exploration 
of the QCD phase diagram. 

The most efficient way to address fluctuations of a system created in a heavy-ion 
collision is via the study of event-by-event (E-by-E) fluctuations, where a given observable 
is measured on an event-by-event basis and the fluctuations are studied over the ensemble 
of the events. Equivalently, one may also study appropriate multi-particle correlations 
within the same region of acceptance [1]. 

In the framework of statistical physics, which appears to describe the bulk properties 
of heavy-ion collisions up to RHIC energies [5], fluctuations measure the susceptibilities of 
the system. These susceptibilities also determine the response of the system to external 
forces. For example, by measuring fluctuations of the net electric charge in a given rapidity 
interval, one obtains information on how this (sub) system would respond to applying an 
external (static) electric fleld. In other words, by measuring fluctuations one gains access 
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Figure 1: Schematic QCD phase diagram. 
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to the same fundamental properties of the system as "table top" experiments dealing 
with macroscopic probes. In the later case, of course, fluctuation measurements would be 
impossible. 

In addition, the study of fluctuations may reveal information beyond its thermody- 
namic properties. If the system expands, fluctuations may be frozen in early and thus 
tell us about the properties of the system prior to its thermal freeze out. A well known 
example is the fluctuations in the cosmic microwave background radiation. Similarly, in 
the context of heavy-ion collisions, one finds that fluctuations of conserved charges of a 
subsystem exhibit relaxation times which are much longer than the typical equilibration 
time of the system [HI E]. 

While the field of event-by-event fluctuations is comparatively new to heavy-ion physics, 
quite a number of measurements have been carried out. Most are concerned with multi- 
plicity fluctuations, fluctuations of particle ratios, transverse momentum, and net-electric 
charge. In this review, we will give an overview of the underlying concepts motivating 
these measurements. While discussing these fluctuations we will keep in close contact to 
results from Lattice QCD calculations, where fluctuations and correlations can be calcu- 
lated for a thermal system. We will also discuss how the careful study of fluctuations may 
help to unravel the QCD phase diagram. In particular, the question about the existence of 
a critical point and its experimental identification will have to involve a detailed analysis 
of the fluctuations and correlations of the system. 

This review is organized as follows. In the following section we will introduce the 
basic definitions and concepts in the framework of a thermal system. Next we will dis- 
cuss fluctuations and correlations of conserved charges. These have a direct mapping on 
Lattice QCD and have the advantage that their relaxation time is rather long, and thus 
should be less affected by hadronic re-interactions. We will then turn to the discussion 
of observables which have actually been measured and try to put them into context with 
recent lattice calculations. The last section will be devoted to a discussion on how the 
study of fluctuations may help in the search for the QCD critical point and the first order 
phase co-existence region. 

Finally, let us emphasize that this article is not intended to be a comprehensive review 
of all aspects concerning fluctuations and correlation in heavy-ion physics. Our main em- 
phasis is on the properties of the bulk matter created in these collisions, which we discuss 
based on a subset of phenomena and observables. For example, multiplicity fluctuations 
will be touched upon only briefly. Also, the interesting correlation structure observed 
among high momentum particles will not be discussed. 

2 Fluctuations and Correlations in a thermal system 
and from Lattice QCD 

The concepts of fluctuations and correlations have a well defined physical interpretation 
for a system in thermal equilibrium. In this case fluctuations and correlations are related 
to the second cumulants of the partition function. These cumulants, or susceptibilities, 
can also be expressed in terms of integrals of equal-time correlation functions, which in 
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turn characterize the (space-hke) static responses of the system. In the case of heavy-ion 
coUisions other, non-statistical, effects may contribute to the measured correlations. For 
example, the dynamical evolution of the system may be too fast for long range correlations 
to build up. These effects need to be understood well in order to extract the matter 
properties from these experiments and we will briefly discuss some of the issues in section 

il 

As already mentioned, the study of fluctuations is essential for the characterization 
of phase transitions. In case of a second order phase transition, the fluctuations of the 
order parameter diverge with a critical exponent specific to the universality class of the 
transition [8j. Furthermore, the system size dependences of the fluctuations can be used 
to distinguish between cross-over transitions and first or second order phase transitions 
using finite size scaling arguments. 

A strongly interacting system in thermal equilibrium can be studied in the framework 
of Lattice QCD (LQCD) (see e.g. [9j for a review). Here we will just use the results from 
LQCD as input to discuss the relevant physic^. 

A system in thermal equilibrium (for a grand-canonical ensemble) is characterized by 
its partition function 



Tr 



exp ( 



where H is the Hamiltonian of the system, and Qi and /Xj denote the conserved charges and 
the corresponding chemical potentials, respectively. In case of three flavor QCD these are 
strangeness, baryon-number, and electric charge, or, equivalently, the three quark flavors 
up, down, and strange. The mean and the (co)-variances are then expressed in terms of 
derivatives of the partition function with respect to the appropriate chemical potential^, 

{Q^) = T^\og{Z) (2) 
{SQ,5Qj) = T^-^-^\og{Z)^VTxi, (3) 

with 5Qi = Qi — {Qi)- Here we have introduced the susceptibilities 

T 

V djiidjji 



XM = T7Tr^M^) (4) 



which are generally quoted as a measure of the (co)-variances. The diagonal suscepti- 
bilities, Xi,ii ^I's a measure for the fluctuations of the system, whereas the off-diagonal 



^We note that at present most lattice calculations are using rather high quark masses, corresponding 
to pion masses of the order of several hundred MeV. However, recently first results with physical strange 
quark mass and light quark masses corresponding to a pion mass of about 220 MeV have been reported 

[H [III nails]. 

^Although in this article we will mostly concentrate on susceptibilities involving conserved charges, 
we note that one can define susceptibilities involving any well defined operator. One prominent example 
is the chiral susceptibility Xch = T /V / dm^^ \og{Z) which characterizes the chiral phase transition in 
QCD. 
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susceptibilities, Xij^ h characterize the correlations between the conserved charges Qi 
and Qj. We note that the susceptibilities are directly related to the well known cumulants 
in statistics [H]. 

One can define and study higher order susceptibilities or cumulants, by differentiating 
multiple times with respect to the appropriate chemical potentials 

1 B'^i B'^j f)'^k 

K,n„n,) ^ ^ ^ J ^ /^N 

^ VTd{fxjTr^d{fx,/Tpd{fXk/Tr^ ^ 

Higher order cumulants up to the sixth [121 IE] even eighth [13] order have been 
calculated in Lattice QCD which, as we will discuss in section [3l provide useful information 
about the properties of the matter above the critical temperature. In Fig. [2] we show the 
Taylor expansion coefficients obtained in two fiavor LQCD [L_5\ which are proportional to 
the susceptibilities defined in Eq. |5l The upper row in Fig. [2] shows the fiavor-diagonal 
susceptibilities up to sixth-order. 



^2 r>T-i9 Xu,U 
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whereas the lower row shows the fiavor off-diagonal susceptibility 

U,d 1 /y\ 

^2 ~ '^^Xu,d {i ) 

and its derivative with the quark number chemical potential fiq = fiu + f^d 

u,d ^ 1 

24T2 9(/i,/T)2^"''' 

U,d _ 1 /ON 

^6 - i44T2 9(/i,/T)4^"'''' ^ ' 

We should point out that the results shown in Fig. [2] are based on simulation with 
rather large quark masses. Recently, new results for three flavor QCD with almost physical 
light quark masses have been reported [10l|16]. In this case, the second order susceptibil- 
ities are consistent with the Stefan-Boltzmann limit of free, uncorrelated massless quarks 
right above the transition temperature T^, while the results shown here (upper left panel) 
exhibit a 20% suppression. The phenomenological consequences of this and other physics 
interpretations of these susceptibilities will be discussed in the following section. 

As already mentioned, susceptibilities are related to intergrals of equal time correlation 
functions of the appropriate charge-densities. Here we will restrict ourselves to the second 
order susceptibilities keeping in mind that the higher order susceptibilities can also be 
expressed in terms of appropriate (higher order) correlation functions. 
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Figure 2: The Taylor expansion coefficients c™ (upper row) of diagonal and c^*^ (lower 
row) of non-diagonal susceptibilities from p[5j. Their relation with the susceptibilities is 
explained in the text, Eq. [6l 



Consider a density fluctuation Spi{x) = Pi{x) — pi at location x, where pi denotes the 
spatially averaged density of the charge Qi. Then the susceptibilities are given by the 
following integral over the density-density correlation functions: 

Xi,j = yf J d^'xd^y {5p,{x) 5pj{y)) = ^pAj + ^ j dhC,,j{r). (9) 

The correlation functions 

c,,(f) = (5p,(f) 5p,(o)) - pA,m - "^""p "^'^'^''^"^ (10) 

are characterized by typical correlation lengths ^jj. The correlation length provides a 
measure for the strength and type of the correlation. For example, in case of a second order 
phase transition the correlation length diverges with a characteristic critical exponent, 
usually denoted as v. 

To illustrate this point let us ffist consider the case of a classical ideal gas. This will 
also serve as useful reference for comparison with LQCD results. Since a classical ideal 
gas has no correlations, by construction its correlation functions vanish, Cideaigas = 0, and 
the susceptibilities are given by the ffist, local term in Eq. |9l ~ PiSi,j, implying that all 
co-variances vanish. As a consequence, the fluctuations are proportional to the number 
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of particles in tlie system, and tlius grow linearly with the system size, V. 

im?) ~ V (11) 

The more relevant case concerning the QCD critical point corresponds to a second 
order phase transition. In this case, the correlation length diverges at the critical temper- 
ature 

^^IT-T^l'" (12) 

where z/ > is relevant critical exponents characterizing a second order phase transition 
in a given universality class |8j. In this case, the volume dependence of the susceptibilities 
is governed by the integral of the correlation function 

x.^^ ~ ~ v'/^ (13) 

so that the fluctuations grow 

{{SQiy) V^/\ secondorder. (14) 

In case of a first order transition we have co-existence of phases with different den- 
sities, and the correlation function is a constant, C{r) = const ^ 0. Consequently, the 
fluctuations scale like 

{Wi?)-V'^, first order. (15) 

Most other systems, including systems with a cross over, such as QCD at vanishing 
chemical potential [17], will exhibit a finite correlation length. Consequently, the suscep- 
tibility is independent of the volume, and the fluctuations scale linearly with the volume, 
just as in the case of an ideal gas 

{(SQiY) ^ V, no phase — transition. (16) 

In principle, one could utilize the above volume scaling of the fluctuations in heavy-ion 
experiments by studying the system size dependence of, e.g., baryon-number fluctuations. 
However, in case of the second order phase-transition, the phenomenon of critical slowing 
down limits the actual size of correlation length due to the finite life-time of the system 
created in these collisions. A maximum correlation length of ~ 2.5 fm has been estimated 
in ref. [18] which is much smaller than the typical size of a system created in these 
reactions. Consequently, such a system would just behave like any other with a constant 
correlation-length and, therefore, would not exhibit the system size dependence discussed 
above. 

After having introduced the basic concepts and definitions for fluctuations in a thermal 
system, in the following section we will discuss, in a few examples, what Lattice QCD 
results on the various susceptibilities can tell us about the structure of the matter above 
the transition temperature. These examples will also guide us in developing observables 
which can and have been measured in heavy-ion experiments. 

■^The correct scaling of the susceptibility with the volume actually involves the critical exponents, 
^ ^ y7/(3i')^ Qm- example here is correct for so called mean field exponents. For details see [8]. 
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3 Fluctuations and correlations of conserved charges 



After having introduced the basic concepts and definitions for fiuctuations and correlations 
in a thermal system, in this section we will discuss some of the physical implications of 
the results obtained from Lattice QCD. Before we do that, however, let us start with 
some general considerations. First, let us address the seemingly contradictory notion of 
fiuctuations of conserved charges. Obviously, if we look at the entire system, none of the 
conserved charges will fiuctuate. However, by studying a sufficiently small subsystem, the 
fiuctuations of conserved quantities become meaningful. The small system may exchange 
conserved quanta with the rest of the system. This is similar to the assumptions which 
govern a thermal system in the grand-canonical ensemble, and Lattice QCD calculations 
are carried out in this ensemblcl. 

dN/dY 




Y 



Figure 3: The various rapidity scales relevant for charge fiuctuations. 

On the other hand, the conservation of the charge ensures that the total charge within 
the small system can only change by transport of charges through the boundaries leading 
to rather long equilibration times for charge fiuctuations in the hadronic phase [H [191 • 
To illustrate this point in the context of heavy-ion collisions, consider a situation as 
depicted in FigJ31 The total system corresponds to all particles distributed in rapidity Y 
over a range Ytotai, whereas the small subsystem corresponds to the particles within the 
accepted rapidity interval Yaccept- For fiuctuations of conserved charges to be a meaningful 
observable the following scales need to be well separated: 

• The range AYfotai for the total charge multiplicity distribution. 

• The interval AYaccept for the accepted charged particles. 

"'The treatment and correction for total charge conservation is non-trivial and needs to be done care- 
fully ^19j since the variance in a grand-canonical and canonical ensemble differ even for large systems 

m- 
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• The charge correlation length AYcorr characteristic to the physics of interest. 

• The typical rapidity shift AY^jcfc charges receive during and after hadronization. 

Given these scales, charge fluctuations will be able to tell us about the properties of the 
early stage of the system, the QGP, if the following criteria are met: 



The first criterion, Eq. [T71 is necessary in order to be sensitive to the relevant physics, 
whereas the second one, Eq. [181 ensures that total charge conservation does not suppress 
the signal and that the signal survives hadronization and the hadronic phase. In particular 
the condition AYaccept ^ AY^ick is unique to conserved charges. For the charge of the 
system to change, charges need to be transported through the boundaries of the system. 
And if the condition AYaccept S> AY^ick is satisfied it requires many kicks to change the 
charge of the system. Consequently the relaxation time into a new equilibrium state may 
be very long, depending on how well the scales are separated. If the charges would not 
be conserved, on the other hand, they could be produced anywhere within the system 
leading to a much more rapid equilibration. These considerations have been explored in 
detail in [7] on the basis of some generalized diffusion equations. 

Hence, for sufficiently high bombarding energies and a sufficiently large acceptance, 
one should be able to see "back" into the QGP, provided that hadronization leads to a 
limited distribution of charge kicks, such that one can assign a typical scale AYkick- The 
charge kick due to re-scattering in the hadronic phase can be estimated using transport 
models and one finds [7] AYkick — 2, close to a naive estimate from p decays, which gives 
AYkick — 1-5. An estimate of the charge transport during hadronization is very difficult on 
the other hand. If, for example, the charge kicks are distributed according to a power law, 
the above considerations are invalid, as the charge transport involves the entire system. 
This might very well be the reason for the results of [25] where in a quark molecular 
dynamics model, the charge fluctuations quickly adjust to the hadron gas value during 
hadronization. We note, that in this model string breaking during hadronization is not 
taken into account, allowing for charge transport over large rapidity intervals due to the 
linearly increasing potential among the quarks. 

Conversely, at low bombarding energies the necessary separation of scales is not nec- 
essarily guaranteed, and the interpretation of fluctuations and correlations of conserved 
charges becomes rather difficult, as charge conservation effects dominate the observables 
[21]. For example, at top SPS energies, Eseam = 160AGeV the width in rapidity of the 
total multiplicity distribution is AYtotai — 2 which is comparable to the charge "kick" 
AYkick > 1.5 due to the decay of the p and Delta resonance. 

After these general considerations about the condition under which the study of charge 
fluctuations in heavy-ion collisions are meaningful let us discuss some speciflc examples 
which have been discussed in the literature. 

The fluctuations of conserved charges in heavy-ion collisions was flrst discused in the 
context of (electric) net-charge fluctuations [ni[2S]- Here the simple observation was that 




(17) 
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charge fluctuations per entropy should scale with the square of the electric charge of the 
charge carrying particles and, consequently, they should be sensitive to the fractional 
charges of the quarks. For simplicity, consider a classical ideal gas of positively and 
negatively charged particles of charge ±g. The fluctuations of the total charge contained 
in a subsystem of particles is then given by 

= q'[{SNl) + {SN^_)] 

= [{N^) + (N-)] = q' (N^h) . (19) 

smce {6N_6N+) = and {6N^) = {N) for an ideal ga^. For a two flavor Quark Gluon 
Plasma this translates into 

(SQl) = Ql (Nu) + Ql {N,) = I {N,) (20) 

where Ng = Nu = Nd denotes the number of quarks and anti-quarks. For a pion gas, on 
the other hand, we have 

{6Ql) = (K.) + . (21) 

Dividing both by the entropy of massless classical particles 

S = 4{N) (22) 

weighted by the appropriate number of degrees of freedom, three in case of the pion gas 
and 40 in case of the QGP, we get 

- 1 (23) 



Sqgp 24 

for a 2-flavor Quark Gluon Plasma whereas for a pion gas we obtain 



(24) 



Thus, the charge fluctuations per degree of freedom in a Quark Gluon Plasma are a 
factor of four smaller than those in a pion gas. Hadronic resonances, which constitute a 
considerable fraction of a hadron gas, reduce the result for the pion gas by about 30 % 
[6[ [27], leaving still a factor 3 signal for the existence of the Quark Gluon Plasma. We 
note that similar arguments can be made for baryon number fluctuations [25], where one 
is sensitive to the fractional baryon number of the quarks. However, in this case an actual 
measurement would require the detection of neutrons, which is rather difficult. 

As already discussed in the previous section both the charge fluctuations and entropy 
can be calculated using Lattice QCD. For the two flavor calculations with rather large 



'"^Corrections due to quantum statistics are small, of the order of one percent for the temperatures 
under consideration [26l [19] . 
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quark masses [15], [9], above the transition temperature both the charge fluctuation and the 
entropy deviate by about 20% from that of a non-interacting gas of quarks and gluons, so 
that the ratio agrees with our estimates. In case of the more recent three- flavor calculation, 
the charge fluctuations are very close to the free gas limit and the entropy deviates by 
at most 10% [llj, so that the ratio again is in good agreement with our simple estimate. 

Unfortunately, at present lattice results are not available for this ratio much below 
the critical temperature. Here one has to resort to hadronic model calculations. This has 
been done in [28l [29] using either a virial expansion, a chiral low energy expansion or an 
explicit diagrammatic calculation. In all cases, interactions slightly increase the ratio as 
compared to a system of non-interacting hadrons, thus enhancing the signal for the Quark 
Gluon Plasma. 

Finally we should note, that present measurements at RHIC and SPS have not found 
the small net-charge fluctuations predicted for a QGP [HI |25] and are rather consistent 
with a resonance gas. 

While the charge fluctuations involve the diagonal susceptibilities, let us next discuss 
an example where the off-diagonal susceptibilities play an essential role. In ref. [3U] the 
following simple observation was made: If above the transition temperature the quarks can 
be well described by uncorrected quasiparticles, then any object which carries strangeness 
necessarily is a quark, and thus carries baryon number. Contrast this with a hadron gas, 
where kaons are mesons which carry strangeness. Thus one expects baryon number and 
strangeness to be stronger correlated in a Quark Gluon Plasma than in a hadron gas. To 
quantify this, let us introduce the following ratio [30] 



C,, . -3-? ^ -3 <^e>-<^,f ^ -3 <||. (25) 



a 



s 



In terms of quark flavors the correlation coefficient Cbs can be written as 

^ o {BS) {{u + d + s){s)) {us) + {ds) Xus + Xds .^^x 

Cb. = -3^ = j^^ = l + + (26) 

since the baryon number of a quark is 1/3 and the strangeness of a strange quark is 
minus one. Note the the quark operators u, d, s here represent the net-quarknumber of 
a given flavor, i.e. {u) = {u — u) etc. Also, to obtain the last expression in Eq. [201 "we 
used the definition for the susceptibilities, Eq. [H from the previous section. For a simple 
quark-gluon plasma, or more generally, for uncorrelated quark flavors, we have 

{us) = {ds) = (27) 

and hence 

Cbs = 1. (28) 

In contrast, a gas of uncorrelated hadron resonances gives 

r ^ o (A) + (A) + ...+3(l]-) + 3(n+) 

- (A^o) + <;^o) + . . . + 9 {Q-) + 9 (0+) • ^^^^ 
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Figure 4: Left panel: The correlation coefficient Cbs for a hadron gas along the empirical 
chemical freeze out line as determined in [5j (full line). The dashed line corresponds to 
and ideal quark-gluon plasma. The figure is from [301. Right panel: Cbs (squares) and 
Cqs (circles) as a function of temperature from Lattice QCD [12j . 

Here the numerator receives contributions only from (strange) baryons, while the denom- 
inator receives contributions also from (strange) mesons. As a result, Cbs = 0.66 for 
T = 170 MeV and /is = 0. 

On the other hand, at very high ^b and low T, where strangeness is carried exclusively 
by Lambdas and Kaons, Cbs ~ f, since (A) = {K) due to strangeness neutrality. This 
significant dependence of Cbs on the hadronic environment is in sharp contrast to the 
simple quark-gluon plasma where the correlation coefficient remains strictly one at all 
temperatures and chemical potentials. In the left panel of Fig. H] we show the result for 
Cbs for both an ideal quark-gluon plasma as well as a hadron gas along the empirical 
chemical freeze-out line, as determined e.g. in [5]. 

The correlation coefficient Cbs has been evaluated in Lattice QCD [12j. The result 
is shown in the right panel of Fig. H] (squares). We see that right above the transition 
temperature Cbs = ^ indicating a system of uncorrelated (quasi)-quarks. Also shown in 
this Figure is the correlation between strangeness and the electric charges 

(QS) {{2u-d-s){-s)) 2xus - Xds 

which also approaches a value of Cqs = 1 right above the transition temperature Tc, 
again indicating uncorrelated quark flavor^. Both results, Cbs — ^ and Cqs — 1 above 
Tc are equivalent, as both result from the fact that the flavor off-diagonal susceptibilities 
are much smaller than the flavor diagonal ones 

Xus < Xss (31) 

Xds < Xss- (32) 



very recent lattice calculation [16j using almost physical quark masses shows a somewhat softer 
transition, reaching the limit of uncorrelated quarks at T/Tc ~ 1.3 
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Figure 5: Negative ratio of off-diagonal to diagonal susceptibility as a function of temper- 
ature in 2-flavor unquenched Lattice QCD. Lattice results are from [15]. Figure adapted 
from 1321. 



The correlation coefficient Cqs may be the better suitable for experimental investigation, 
since, contrary to Cbsi it does not involve the measurement of neutrons. There may be 
even better combinations of quantum numbers as discussed in [31] . However they all rely 
on the fact that the ratio of flavor-off diagonal to flavor diagonal susceptibilities is small 
above the transition temperature. 

A similar argument can be made using only light flavors. In this case one would look 
at something like baryon-number - "up-ness" correlations. Again the relevant quantity 
would be the ratio of flavor off-diagonal to flavor diagonal susceptibilities, Xud/Xdd and 
uncorrelated flavors would correspond to 



Xud 
Xdd 



< 1. 



(33) 



In FiglSl we show this ratio for the light quarks based on the Lattice data of ref . [15] . 
The diagonal and off-diagonal susceptibilities are also shown separately in Figj2l We find 
that the ratio drops rapidly as we approach the transition temperature T/Tc = 1 and for 
temperatures T /T^ > 1.2, the ratio is consistent with zero, indicating independent quark 
fiavors. This behavior is consistent with a rapid melting of the hadronic states in the 
transition region. It also suggests that the system above T^. should be well described in 
terms of independent quasi particles. And indeed, many of the lattice results, including 
the susceptibilities can be reproduced in quasiparticle [53] and mean field models [5^ . 

Further evidence for independent quasi-quarks is found in the higher susceptibilities. 
In [35] the authors looked at the ratio of the fourth order susceptibilities (cumulants) over 
the second order susceptibilities for baryon number and net charge 



R 



iB,Q) 



'4,2 



(4) 
(2) 



(34) 



where we used the notation of EqlS] In Figl6] we show lattice results from a recent 
calculation with three fiavors and almost physical light quark masses [TOl [16]. On the 



left panel we have the ratio i?f 2 baryons and on the right panel the ratio R'^2 the 
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electric charge. In both figures the full lines indicate estimates for a hadron gas (labeled 
"HRG" ) at low temperatures and the limit of non-interacting quarks at high temperatures, 
(labeled "SB). In addition, the left panel shows two estimates for the hadron gas, one with 
(higher) and one without pions (lower). We note that for both ratios the peak close to 
the transition temperature softens considerably when going to smaller lattice spacings, 
i.e, from A^^ = 4 to A^,- = 6. 

Let us now concentrate on the left panel, where the ratio for the baryon number 
susceptibilities is shown. This result can be easily understood in terms of hadrons on the 
low temperature side and independent quarks on the high temperature side. Consider a 
classical ideal gas of particles with baryon number b. Then we have 

xf = b-(x'^ + xf) (35) 
Xf = f'(x!J'+X^') (36) 

(2) (2) 

where Xn ^^"^ x)^ s-^'e the particle-number cumulants for particle and anti-particles, 
respectively. 

x!?^ = {N')-{Nf = {N) (37) 
xS^ = {N') ~ 3 {N') = (N) (38) 

Consequently 

jB _ l2 



i?f2 = ^ (39) 



and since all baryons in the hadronic phase have baryon number \Bhadronic\ = 1 and all 
quarks have baryon number \Bquark\ = 1/3 we get for the the ratio of the cumulants the 
simple results 

= 1; hadron phase (40) 

= I; QGP (41) 

(42) 

Since baryons are fermions, one would have to correct the above result for quantum 
statistics. In case of massless particles this can be done analytically and one would have 
to multiply the results by a factor of G/vr^ ~ 0.6. The effect of quantum statistics for the 
massive baryons in the hadronic phase can only be evaluated numerically and for baryons 
with mass M = 1 GeV and a temperature of T = 170 MeV the correction is less than 1%. 
The results of this simple estimate are also shown in FigE] as the full lines at low (hadron 
gas) and high (QGP) temperatures. 

Again, we find that above ~ 1.5 Tc the system is well described by a gas of independent 
quarks. Thus, the study of correlations and fiuctuations reveals insights into the structure 
of strongly interacting matter above the transition temperature which the equation of 
state alone could not tell us. In addition, if the latest A^,- = 6 with light quarks are 
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Figure 6: Ratio of fourth order to second susceptibilities for baryon number (left panel) 
and electric charge (right panel). The lines at low temperature, labeled "HRG" indicate 
the results for a hadron gas. The limit of free quarks is shown on the right of each plot 
and denoted by "SB" . The two hadron gas lines on the left panel represent a hadron gas 
calculation with (upper) and without pions (lower). Figure adapted from P^ . 

indeed what the continuum limit will look like, then we see a rapid change from hadrons 
to independent quarks, both in the flavor off-diagonal susceptibilities as well as in the 
fourth order baryon number cumulants. This would lend theoretical support of the rather 
surprising finding of the so-called quark number scaling at RHIC [36], [371 EH] , which can be 
simply understood within a recombination/coalescence picture. In this picture hadrons 
are formed at by simple phasespace coalescence, which is consistent with the the rather 
rapid buildup of correlations around seen in the lattice results for the susceptibilities 
as discussed here. If, on the other hand, the peak in the fourth order baryon-number 
cumulant just below of the A^,- = 4 results is what the continuum limit will be, then 
yet to be understood strong correlations are at work close to the transition temperature. 
This in turn could be a sign of the existence of a critical point at finite baryon-number 
chemical potential, where all the baryon-number cumulants are expected to diverge. The 
physics of the critical point will discussed in more detail in section [51 

To summarize, in this section we have shown that the study of fluctuations and corre- 
lations reveals important insights into the underlying structure of the strongly interacting 
matter, especially above the transition temperature. How this can be translated into 
experimentally accessible observables will be the subject of the following section. 

4 Observables 

In the previous sections we have concentrated on fluctuations and correlations of thermal 
systems and have presented a few examples on how to extract interesting physics from 
them. Now let us turn to actual heavy-ion collisions and let us discuss how to construct 
meaningful observables involving the various susceptibilities discussed previously. 

One essential difference between a system studied in lattice QCD and a heavy-ion 
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collision is that in the latter case the observables are restricted to correlations in the 
momenta and quantum numbers of the observed particles. Spatial correlations, such as 
the correlations length, are accessible only indirectly via Fourier transforms of momentum 
space correlations, and thus limited. The susceptibilities, on the other hand, can be 
measured since they can be expressed as integrals over either spatial or momentum space 
correlation functions. Thus, as long as one deals with susceptibilities, i.e (co) -variances, 
there is a one to one mapping from Lattice QCD results to heavy-ion collisions (assuming 
that the latter produce a system in thermal equilibrium). The susceptibilities can be 
extracted from data either by studying event-by-even fluctuations of a given quantity or 
by measuring and integrating the appropriate multi-particle densities [4]. 

An additional complication in the experimental study of fluctuations in heavy-ion 
collisions are unavoidable impact parameter fluctuations. These additional sources for 
fluctuations may very well mask the fluctuations of interest. In the language of thermal 
systems, these impact parameter fluctuations correspond to volume fluctuations. Conse- 
quently, one should study so called intensive variables, i.e. variables which do not scale 
with the volume, such as temperature, energy density, etc. 

To simplify the discussion, let us concentrate on two-particle correlations, which fully 
characterize fluctuations of Gaussian natur^. Let us begin by defining the number of 
particles with quantum numbers a = {«!, . . . , in the momentum space interval {p,p + 
dp) in a given even@ 

dN" 

. = _pst dp (43) 
^ dp 

and its fluctuations 

6n; = n^- (n;) (44) 
where (^p) denotes the event-averaged multiplicity in the momentum bin {p,p + dp) 

event j=l 

In the case of an ideal gas, (^p ) takes the form of the Bose-Einstein or the Fermi-Dirac 
distribution depending on the spin of the particles. 
The mean value of an observable 



X = Y.^'^n'; (46) 



p,a 



^In the case of fluctuations of positive definite quantities, such as transverse momentum or energy, 
the appropriate distribution is a so called Gamma- distribution [39] 

^In this section, we will use the term 'event' and 'member of the given ensemble' interchangeably. We 
also use 'event average' and 'ensemble average' interchangeably. 

^Here we use sums over momentum states, as appropriate for a finite box. The conversion to continuum 
state is a usual, J2p ^ ^ I (t-J)'^ ' where V is the volume of the system. 
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is obtained by averaging over all the events in the ensemble 



p,a 



Note that X defined in this way is an extensive observable and thus subject to volume 
fiuctuations. 

For the discussion of fiuctuations and correlations we need to consider two-particle 
densities. Let us, therefore, define the basic correlator A,"'f 



p,<} 



Ap,f = (Sn^M) (48) 

where 6n'^ is defined in Eq.f l44p . In case of an ideal gas, particles with different quantum 
numbers and in different momentum states are uncorrelated. Therefore, 

K:S = ^p'^ (fl^-^p^ ^pi^p) (49) 

with 

< = 1 (50) 



for a classical ideal gas, and 



a 



UJp 



:i ± K» (51) 



for a boson (+) and a fermion (— ) gas (see e.g [iQj). In the presence of dynamical corre- 
lations the basic correlator A^'^, EqH8| may contain off-diagonal elements in momentum 
space and/or in the space of the quantum numbers. For example, the presence of a res- 
onance, such as a po, not only correlates the number of vr"*" and vr" in the final state due 
to resonance decay, but also their momenta. In addition, the occupation number Up in a 
given momentum interval may be changed, for example, due to collective fiow effects. 

Any two-particle observable can be expressed in terms of the basic correlator A^'^ 
[H m] . For the generic observable X as given by Eq. SHI the variance is given by 

{6X') = J2 ^pi« (52) 

To illustrate the formalism let us consider the net charge 

g = 5^g(«)n^. (53) 

Here, Q{a) is the charge of the particle with the quantum numbers a = {ai, . . . , Its 
variance, i.e. the net charge fluctuations, are given by (see Eq{ 



[SQ') = Yl Q{o^)Q{(^)Ki- (54) 
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In case the system contains only particles with (positive or negative) unit charge, the 
above expression simplifies tc0 

P,Q 

= {5Nl) + {5N^_) - 2 {5N+5N^) . (55) 

Note, in the last equation the basic correlator still contains a sum over all other 

quantum numbers, which we have suppressed. 

Kt= E Ai;^-M+A} (56) 

ai T^charge ,/3i T^charge 

and similar for the other combinations, meaning that all charged particles are included 
independent of their fiavor, spin, etc. 

For a thermal system, the basic correlator A^'^ is directly related to the second order 
susceptibilities 

{SQi6Qj) = J2 Q^i»)QM A°f = VTxi,. (57) 

Here, Qi{a) denotes the conserved charge Qi, such as electric charge, baryon number, or 
strangeness, associated with the quantum number a = {ai, . . . For instance, the 

baryon-number strangeness correlations discussed in section [3] are given in terms of the 
basic correlator as 

{BS)= J2 Bia)S{P)A;f. (58) 

Obviously, all the information revealed by (Gaussian) fiuctuations is encoded in the basic 
correlator A^'^ and the fiuctuations of different observables expose different "moments" 
and elements of the basic correlator. As we will discuss below, all the event-by-event fluc- 
tuation observables studied so far, such as ratio fluctuations, mean transverse momentum 
fluctuations, etc, can be expressed as moments of the basic correlators. 

If the system has no genuine two-particle correlations, then the basic correlator will 

be 

A^,f = {n,) 5,,,S^,p (59) 

similar to the classical ideal gas, except that (n^) does not have to follow a Boltzmann 
distribution. Using the inclusive single particle spectrum as momentum distribution, the 
relation Eql59] may be utilized as an estimator of the fluctuations due to flnite number 
(Poisson) statistics [52]. A more detailed discussion can be found in [T9] . 

Any correlations, on the other hand will lead to a non Poisson behavior of the basic 
correlator. The ideal quantum gases discussed above are examples where the symmetry 

^"^Note the difference in notation: whereas Snp refers to the fluctuations in the momentum interval 
{p,p + dp), SN refers to the fluctuations of the total (integrated) number of particles. 
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of the wavefunction introduces 2-particle correlations which either reduce (Fermions) or 
enhance (Bosons) the fluctuations. Also global conservation laws such as charge, energy, 
etc. will affect the fluctuations. 

As mentioned in the beginning of this section, in order to avoid contributions from 
volume / impact parameter fluctuations it is desirable to study so called intensive quan- 
tities, i.e. quantities which do not scale with the size of the system. As already pointed 
out, in a heavy-ion experiment all observables need to be constructed from the momenta 
and quantum numbers of the final state particles in each event. Furthermore, the number 
of produced particles typically increases with the centrality of the collision, i.e. with the 
volume of the system created. Therefore, the only way to devise volume independent ob- 
servables is by considering ratios of observables which have the same volume (centrality) 
dependence. Examples which have been explored in experiments are ratios of particle 
numbers, mean energy or mean transverse momentum, i.e energy or transverse momen- 
tum divided by the number of particles. Thus, it will be useful to discuss the fluctuations 
of ratios in some more detail. 

4.1 Fluctuations of Ratios 

Consider the ratio of of two observables A and B 

Rab = 4 (60) 

witlil^ 



5 = Ew- (61) 



Expanding A and B around the mean, A = (A) + SA and B = (B) + 6B, we get 

«,,.M,i4>fiA_i£Vo(^^) (62) 



(B) (B) \(A} (B> 

where 0(5^) indicates terms that are of quadratic and higher orders in 6 A and 6B. Since 
A and B are extensive observables, the neglected terms are at most 0{1/ (N)). From 
Eq. (p21) . it is easy to see that to leading order in 6 

(Rab) = 1^ (63) 

and the variance 

^ / {Bf V {Af {Bf {A) (B) J ^ ^ 



^^Defined this way, both A and B scale with N, the number of particles in the final state. 
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Using the basic correlator A^'^, this can be rewritten as 

((ifi.«)>-^j:j:A,,^^+^-2^j. (65) 

Often in the hterature {5R\g) is scaled by the square of the mean ratio to obtain a relative 
width, 

which in case of an ideal gas simplifies to 

ai,(,deal gas) = E (°S> < + ^ " ^ (||) ) ■ (67) 

In addition, covariances between two different ratio can be expressed in terms of the 
basic correlator in similar fashion 






{A) (C) 




(B) (D) 







{A){C) {B){D) {A){D) {B){C) 



where 



c = E«' ^ = E«- (69) 



This becomes useful, for example, in the discussion of charge dependent and charge inde- 
pendent transverse momentum fluctuations |43j . 

The interesting physics, of course, resides in the dynamical correlations of the system. 
In order to expose these, one needs to compare to the purely statistical fluctuations of an 
uncorrelated system, which experimentally can be obtained via mixed event techniques. 
For ratio fluctuations, the variance for an uncorrelated system is given by the expression 
of the ideal gas, Eql67]with a;^ = 1, and (n^) corresponding to the measured inclusive 
momentum distributions of the particles with quantum numbers a. We can define the 
uncorrelated variance CFl^r^coTrcUted 



(^uncorrelated - 2^\^p)^p ^ ,^>2 + ,^v2 ^ ^^^^ j 



(70) 



P 

where here (...) corresponds to actually measured (inclusive) averages. The difference 
between the actual variance of the system and that of an uncorrelated (mixed event) 

dynamic 



ensemble is usually referred to as 



^ dynamic ^ ^ uncorrelated " V / 
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The first such measure to be proposed has been the so called $ variable |l2l |15] which in 
terms of our variables here is given by 



= (Rab) ^/{N){a - (Tuncorrelated). (72) 

The ratio has also been proposed 



a 



(73) 



uncorrelated 



which has the advantage of being dimensionless and being less sensitive to acceptance 
effects, as we shall discuss in some more detail in the context of the fluctuations of the 
kaon-to-pion ratio below. We finally note, that with these conventions, correlations due 
to quantum statistics already contribute to e.g. adynamic; typically at a few percent level, 
as discussed in HH Il5] . 



4.1.1 Examples: Fluctuations of the mean transverse momentum and parti- 
cle ratios 

Event-by-event fluctuations of the mean transverse momentum, pt, have been discussed 
in the literature as a measure of energy fluctuations, which should show a peak close to 
the QCD phase transition, where the specific heat has a maximum [l6l HTj. Transverse 
momentum fluctuations have also been discussed in the context of a search for the QCD 
critical point. As we will discuss in more detail in section [5l close to the critical point 
one expects long range fluctuations which would result in enhanced transverse momentum 
fluctuations, especially for small momenta [HI HH] . The signature in this case would be a 
maximum in the excitation function of pt-fluctuations at the energy corresponding to the 
location of the critical point. Before we discuss some of the experimental results, let us 
first write down the necessary definitions. The mean transverse momentum is defined as 

where (Pj) denotes the transverse momentum of an event with particles in the final 
state. Obviously, we are dealing with a ratio and with the substitution 

A = Pt, ttp = sj pI + pI = Ptip) 

B = N; b^ = l (75) 

the above formalism, Eqj66l can be readily applied leading to the expression for the scaled 
variance 

2 ^a,l3 f Pt{p)Pt{q) , 1 o Pt{p) \ (^p_. 
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which in the case of an ideal gas (see Eqj671) reduces to 



(ideal gas) = ^ (n;^) (^pt(p) - j 



(77) 



We note that, since (Pt) scales with the number of particles in the final state, {Pt) ~ (A^), 
the variance scales as a^^ ~ j^, and thus is sensitive to the actual acceptance of the 
detector. The same is true for all ratio fluctuations we discuss here. 

Experimentally, pt-fiuctuations have been investigated by the CERES collaboration 
p], the NA49 collaboration [SDl El] at the CERN SPS and by the STAR |13], and the 
PHENIX [52| [53] collaborations at RHIC. These measurements cover a wide range of 
beam energies. The resulting p^-fluctuations are shown in Figd They are small and show 
virtually no beam-energy dependence. 
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Figure 7: Transverse momentum fluctuations for different center of mass energies. Here 
Spy = (Tdynamicai [5^- The figurc is adapted from ^54j. 

In addition to the transverse momentum fluctuations for all charged particles, one can 
investigate the pt fluctuations of the negative and positive charges independently as well 
as the cross correlation between them. This aspect is discussed in detail in \19\ and has 
been investigated in experiment [l3l[55]. In addition to the charge dependence the NA49 
collaboration has also analyzed the pt fluctuations for different cuts in the transverse 
momentum (see Fig JTU]) in order to test ideas about the QCD critical point, which we will 
discuss in sectional 

Let us next turn to the fluctuations of particle ratios and specifically to the fluctuations 
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of the kaon-to-pion ratio, K/n. In this case 



B 



A 




(78) 



so that 




2 



WW 



) 



(79) 



In the absence of any correlations, i.e a classical ideal gas, this reduces to 



c^i/^(uncorrelated) = — + — 




(80) 



The measurement of the K/tc fluctuations by the NA49 collaboration [56] were the first 
event-by-event fluctuations measurement in a heavy- ion experiment. The original mo- 
tivation for this analysis was to look for separate event classes, for example, one with 
enhanced strangeness and one without. The observed fluctuations, however, are rather 
small, ~ 2%, indicating that the events generated in these collision are rather similar. 
Subsequently, the NA49 collaboration has measured the the K/n fluctuations over the 
entire range of energies available at the CERN SPS |57], and together with the prelim- 
inary measurements from STAR at RHIC [SS] we have an excitation function for this 
observable over a wide range of energies. This is depicted in the upper panel of FigJHl 
We observe a rather steep increase of the K/n fluctuations as the center of mass energy 
is decreased below \fs ~ lOAGeV. This rise coincides with a maximum of the inclusive 
K/tt ratio [59]. It can neither be explained in the statistical hadron gas model [60] nor 
can it be reproduced by the URQMD model [61] (see FigJHD, which is able to describe the 
drop in the fluctuations of the proton-to-pion ratio, shown in the lower panel of Fig. [HI 

Consequently, this rapid increase has sparked considerable discussion and speculations 
concerning the QCD critical point and the first order co-existence region, which we discuss 
in more detail in section [51 While this might indeed be a first experimental indication 
for interesting structures in the QCD phase diagram, we would like to raise a note of 
caution. As discussed in [27] , ratio fluctuations scale roughly as the inverse of the accepted 
multiplicity. 



Consequently, the observed rise may partially be due to the change of the actual accep- 
tance with beam energy, which is always the case in a fixed target experiment such as 
NA49. As proposed in [27], the observed rise in the fluctuations would be more convinc- 
ing if it showed up in the ratio of the measured variance over that from mixed event, as 
defined in Eq. [731 This ratio does not exhibit any trivial multiplicity dependence. 

Finally let us briefly discuss the fluctuations of the ratio of positively over negatively 
charged particles, which is directly related to the net-charge fluctuations discussed in 
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Figure 8: Upper panel: Fluctuations of the K/n ratio as a function of the center of mass 
energy. The data at low -^i are from the NA49 collaboration [57j and the ones at high -^i 
are prelimnary data from the STAR collaboration [58j. Lower panel: Fluctuations of the 
p/tt ratio as a function of the center of mass energy from NA49. Results from URQMD 
|6T] calculations are shown as open squares. The figure is adapted from [57] . 



section El These fluctuations are a direct probe for the existence of a QGP as pointed 
out in [ni 125] and the many subtleties such as global charge conservation [B2] have been 
discussed in the literature (for a review see [19]). The fluctuations of the ratio of positively 
over negatively charged particles 



is given by 



R.- = ^ (82) 



which in the limit of small net charge 

(Q) = {N+ - N_) « {N,h) = (iV+ + iV_) (84) 
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can be written as 




(5Nl + 5N^ - 2SN+SN^ 



4 m 



(85) 



Since the number of charged particles is a measure for the entropy, (Nch) ~ 5*, the 
observable D = (Net) measures the charge fluctuation per degree of freedom, which, 
as discussed in section [31 changes by roughly a factor of three when going from a hadron 
gas to the a Quark Gluon Plasma. As already mentioned, for this observable to be 
useful, a clear separation of scales in rapidity is necessary, otherwise hadronization and 
subsequent hadronic interactions degrade the signal. This observable has been measured 
in experiment both at the SPS and at RHIC |SSl El ESI EH] by the CERES, NA49, 
PHENIX, and STAR collaborations. While the RHIC results are consistent with the 
predictions for a resonance gas, the results at the SPS are difficult to interpret as the 
necessary separation of rapidity scales at these energies is not satisfied and global charge 
conservation effects dominate the signal. No indication of a Quark Gluon Plasma has 
been seen, however. 

5 Fluctuations and the QCD critical point 

As pointed out in the introduction, fluctuations and correlations are unique signatures for 
phase transitions. Therefore, an experimental search for a possible critical point and a 
first order co-existence region in the QCD phase diagram is intimately connected with the 
study and measurement of fluctuations and correlations. In this section we will discuss 
how presently available data constrain the existence of a critical point and what additional 
measurements should be done to shed more light on this very interesting question. 

Let us start with a brief reminder of the physics and the present theoretical under- 
standing of the QCD critical point. A detailed account can be found in [HTj. The basic 
argument which leads to the existence of a critical point in the QCD phase diagram is 
rather straightforward. For vanishing baryon-number chemical potential Lattice QCD 
calculations have established that the transition is an analytic cross-over [17J. If the tran- 
sition at zero temperature but flnite density is of flrst order, as most models predict, then 
the flrst order phase transition line in the T — /iBaryon plane has to end. This "end-point" is 
the QCD critical point associated with a second order transition. This scenario is depicted 
in Fig. [H The question of course arises how robust this simple argument really is. The 
cross-over transition at /is = seems well established by Lattice QCD. The flrst order 
transition at T = 0, on the other hand, is model dependent. While essentially all chiral 
models do predict such a flrst order transition at T = and flnite /ibaryon one should keep 
in mind that all these modelfl are essentially extensions of the linear sigma model, such 
as the Nambu model with or without Polyakov loop dynamics, and small modiflcations 
may already alter the conclusions. For example, once the t'Hooft anomaly interaction is 
taken into account, the flrst order transition at T = may turn into a cross-over with yet 
another critical point at low temperatures and high density |68] . 

^^For a compilation of the model predictions see [67] . 
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A nice way to analyze these models and also Lattice QCD results in the region of 
small chemical potentials is to find the critical quark mass for which one obtains a second 
order transition [69]. This is depicted in Figl9l Most chiral models predict that the 
critical quark mass increases with the chemical potential (right panel of FiglH]). In this 
case, one expects a critical point at finite chemical potential once the critical quark mass 
coincides with the physical quark mass, as can be seen in the figure. Lattice QCD, on the 
other hand, seems to favor the opposite trend, namely a decreasing critical quark mass 
(left panel of Fig 19]). This is the result of [691 HO] obtained in an expansion up to fourth 
order in the chemical potential /i^ on a rather small lattice. If these lattice results hold 
up for larger lattices it seems that the chiral dynamics does not predict the small /i^ 
behavior of the critical quark mass correctly and other effects are more dominant. One 
possibility would be a repulsive vector coupling, which is neither constrained nor ruled 
out by symmetry arguments. As shown in [71] a suitable choice of a repulsive vector 
coupling can indeed reproduce the trend seen on the lattice. 

There have also been attempts to find the critical point directly in Lattice QCD. 
This requires reweighting techniques in order to circumvent the problems associated with 
the complex fermion determinant at finite chemical potential. The pioneering work [72] 
in this approach has indeed located a critical point. With realistic quark masses, this 
method predicts its location at T ~ IGOMeV and Hb — 360 MeV [73]. However, the 
method employed can not easily be extended to larger volumes and, therefore, one does 
not know if the signal survives in the infinite volume limit. Other approaches calculate 
the free energy at finite chemical potential as a Taylor expansion in terms of the chemical 
potential (see e.g [15], [13]). As discussed in section [31 the expansion coefficients are given 
by the baryon number cumulants or susceptibilities. While this method does not allow 
to extract a critical point directly it can provide limits for its location. At present a 
conservative limit for its chemical potential is /x^ > 2Tc ([HI US]), where ~ 180 MeV 
is the transition temperature at vanishing baryon number density. 



Figure 9: Behavior of the critical line as a function of chemical potential hb- Left panel: 
Scenario favored by Lattice QCD [691 170] where critical line moves towards smaller quark 
masses. Right panel: Standard scenario predicted by most chiral models, where critical 
line moves towards higher quark masses. The figure is adapted from [69j. 
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Obviously there is only limited theoretical guidance for an experimental search for the 
critical point as the model predictions for its location vary quite a bit. From hadronic 
freeze out systematics [5J, on the other hand, one knows that the chemical potential of 
the system created depends on the center-of-mass energy of the collision. Unless the 
temperature of the critical point is unexpectedly low, one can explore regions up to 
about /iB < 500 MeV in the chemical potential by lowering the beam energy to about 
^/s ~ 5 GeV. Hence, the strategy for a search is to study excitation functions of vari- 
ous observables and see if they show non-monotonic behavior at the same beam energy, 
indicating the location of the critical point or of the first order phase co-existence region 

Since the baryon density is an order parameter for the phase transition at finite density, 
baryon number fluctuations are the natural observable to consider. In section |3] we have 
already discussed the serious limitations baryon number conservation imposes on this 
observable, especially for low center-of-mass energies. In addition, the measurement of 
the baryon number requires the detection of neutrons, which is difficult. However, as 
argued in [73], it may be sufficient to study proton number fluctuations, as the iso- vector 
channel does not show critical behavior. This may also soften the limitations due to global 
baryon-number conservation. However, even if the system reaches the critical point, the 
correlation length, which diverges in a thermal system, would be finite due to critical 
slowing down together with the finite time the system has to develop the correlations. In 
^18j a correlation length of ^ ~ 2.5 fm has been estimated based on these considerations. 
Therefore, it may be advantageous to study higher order cumulants which depend on 
higher powers of the correlations length [75j. Indeed the fourth order cumulant x^^\ 
Eq. |38|, scales like the seventh power of the correlation length, ~ [75]. Thus if the 
correlation length increases only by 10% in the vicinity of the critical point, one should 
see an enhancement by a factor of two in the fourth order cumulant, whereas the second 
order cumulant, i.e. the fiuctuations, would only increase by 20%. 

Initially transverse momentum fluctuations have been proposed as a signature for the 
critical point [l8lllT|, since close to the critical point the system should develop large, and 
mostly long range, i.e low momentum, fluctuations. Therefore, it was suggested that an 
excitation function of the transverse momentum fluctuations should show non-monotonic 
behavior, especially for small transverse momenta. In the meantime, such an excitation 
function has been measured and it is shown in Fig {10] for different charge combinations and 
different cuts on the transverse momentum. Critical fluctuations, corrected for critical 
slowing down and expansion of the system [T8| would lead to a bump which should be at 
least a factor of two larger than the statistical background. Obviously, the data shown in 
FigHni do not show such a behavior, even for small transverse momenta. The results at 
RHIC [121, shown in FiglTJ are consistent with the data from SPS. Hence, so far there is 
no indication of a critical point in the transverse momentum fluctuation measurements. 
Of course it could be that the signal is too weak to be seen and it may also be washed 
out by subsequent hadronic interactions. To address this issue, higher cumulants, as 
discussed above, need to be measured as they should show a stronger enhancement close 
to the critical point. Furthermore, on the theoretical side, one needs to get a better 
understanding of the degradation of the proposed signals in the hadronic phase. 

The only observable which shows a strong beam energy dependence are the fluctuations 
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Figure 10: Preliminary data on the energy dependence of pt fluctuations from the NA49 
collaboration [55] for all charged particles and for positively and negatively charged par- 
ticles. The panels show the fluctuations for different cuts in the transverse momentum. 
The figure is adapted from [55] . 

of the kaon-to-pion ratio, shown in Fig|Sl As discussed in section HJ neither a transport 
approach nor the statistical hadron gas model can reproduce these data. However, as 
already pointed out, o"^ynamic scales with the inverse accepted multiplicity, and the observed 
rise may very well be partially due to the changing acceptance of the fixed target NA49 
experiment. If the enhancement turns out to be real, it is still not clear if is related to 
the critical point. Close to the critical point one would expect the fluctuations of the 
pion number to be enhanced. But this would imply also enhanced fluctuations of the 
proton-to-pion ratio, which is not observed in experiment, as shown in the lower panel of 
FigEl There the data are well reproduced by the URQMD calculations. 

Although this section is mostly about the critical point and its detection, let us em- 
phasize that it might be more beneficial to look for and identify the first order co-existence 
region. Contrary to the critical point, the first order transition corresponds to an entire 
region in the T — p phase diagram. Thus it is more likely for the system to cross this re- 
gion rather than the critical point. It is also more likely for the system to spend sufficient 
time in this region in order to develop measurable effects. One example is the develop- 
ment of spinodal instabilities, which are a generic phenomenon of dynamical first order 
transitions [75] . Spinodal instabilities have been studied and successfully identified in the 
context of the nuclear liquid gas phase transition [77]. In the case of the QCD first order 
transition, spinodal instabilities could lead to kinematic correlations among particles [_78j 
and to enhanced fluctuations of strangeness [79]. And indeed the observed enhancement 
of the kaon-to-pion fluctuations, if real, may be due to these enhanced fluctuations in the 
strangeness sector [70]. However, just as for the critical point, there has been no quantita- 
tive calculation of the effect due to hadronic re-scattering on the observables. In addition, 
so far there is no dynamical model which carries the system through the spinodal region, 
as it is the case for the nuclear liquid gas phase-transition, where such as model proved to 
be extremely useful in guiding the experimental searches. These models have also helped 
to develop unique observables, such as the variance of the cluster size which subsequently 
could be identified in experiment |77j. This analysis lead to a rather convicing case for the 
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existence of a first order phase co-existence region for nuclear matter. 

Although this article is about fluctuations, in closing let us briefly mention other 
observables which are discussed in the context of the critical point. Most prominently 
is the idea to look for soft modes in the low-mass dilepton invariant mass spectrum. 
However, it is not clear if the soft modes, responsible for the large density fluctuations 
close to the critical point, are visible in the dilepton channel, since they are of space like 
origin. Indeed, an analysis of the fluctuations close to the critical point carried out in 
the Nambu model with finite quark masses [801 E] shows that the sigma-meson remains 
gaped at the critical point, contrary to the chiral transition in the limit of vanishing quark 
masses. Thus, in this model no significant change due to the critical point has been seen 
in the time-like spectrum, which is accessible to dilepton spectroscopy. 

There are a number of other possibilities which have not yet been explored theoreti- 
cally. For example, it maybe interesting to further explore the co-variances between the 
baryon density fiuctuations and other quantities which couple to the baryon density, such 
as e.g. dileptons. These co-variances are expected to become large and could possibly 
be developed into practical observables, which will not be affected by baryon number 
conservation. 

6 Summary and Conclusions 

In this article we have provided a review on some selected aspects concerning the study of 
correlations and fiuctuations in heavy ion collisions. We have discussed some of the recent 
lattice results on fiuctuations and correlations. While the results are still improving, a 
rather interesting picture for the matter above the transition temperature emerges from 
these Lattice calculations. The ratio of fiavor off-diagonal to diagonal susceptibilities 
as well as the fourth order cumulants in baryon number and charge indicate that the 
quark degrees of freedom behave like uncorrelated particles at temperatures above T > 
1.3 Tc. Furthermore, the transition from a hadron gas seem rather rapid. To which extent 
there are some additional, non-trivial features just below the transition temperature is, 
at present, an open question. While the most recent Nf = 6 data with almost physical 
quark masses show a rather featureless transition from hadron to independent quarks, the 
results for a coarser, Nt = 4, lattice exhibit a peak in the fourth order cumulant for the 
baryon number. This, if correct, would be a very nice first indication for a critical point 
at finite chemical potential. If, on the other hand, this peak disappears in the continuum 
limit, as the Nt = 6 data seem to suggest, the critical point will likely be located at much 
higher chemical potential. We note, that the featureless, rapid transition from hadrons 
to quarks, would nicely explain the observed quark number scaling and recombination 
phenomenology observed at RHIC. 

As far as observables are concerned we have discussed electric charge fiuctuations, 
transverse momentum fiuctuations, and the fiuctuations of the kaon-to-pion ratio. All 
these have been measured over a wide range of beam energies, from the CERN SPS to 
RHIC, and none of the excitation functions, with the possible exception of the kaon-to- 
pion ratio, show any significant beam energy dependence. In the case of the kaon-to-pion 
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ratio a rapid rise towards the lowest energies is seen, which may or may not be due to 
simple scaling of the observable with the acceptance. If the observed rise is real, then this 
may very well be the first hint for some non-trivial phase structure probed at these lowest 
energies. 

Concerning the QCD critical point, the present data set does not provide any evidence 
for its existence in the region probed so far. The transverse momentum fluctuations do 
not show any non-monotonic behavior, as originally predicted. However, it could very well 
be that the signal is too weak, as there is not sufficient time to develop a large correlation 
length in these systems. Therefore, it is imperative to measure the higher cumulants 
as well. For example, the fourth order cumulant scales like the seventh power of the 
correlation length, whereas the second order, which controls the transverse momentum 
fluctuations, only scales like the square of the correlation length. Fortunately, the RHIC 
accelerator can easily cover the range of energies required for a comprehensive exploration 
of the QCD phase diagram. And, with the two state of the art detectors, PHENIX and 
STAR, a definitive measurement of these and other observables is possible. In addition the 
CBM experiment at the FAIR facility will be able to perform high statistics measurements 
in the region of high chemical potential, allowing for a detailed spectroscopy of this most 
interesting region of the QCD phase diagram. 

Let us conclude by pointing out the obvious: A dedicated energy scan program at 
RHIC and FAIR would be very desireable to definitively answer the central question of 
the field: Is there a non-trivial phases structure of QCD? 
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